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Theory 

(1)  Electronic  Structure  of  Disordered  Systems. 

The  study  of  the  electronic  structures  of  one-dimensional  disordered 

arrays  of  delta-function  potentials  has  provided  us  with  considerable  insights 

into  the  applicability  of  the  Wu-Dy  method (1)  to  three-dimensional  systems. 

The  preprint  of  a  paper  containing  the  results  is  attached.  Extension  to  three 

dimensions  are  now  under  way.  The  models  considered  are  those  of  Henderson(2) 

(3) 

and  Polk  .  The  function  R^,  required  in  the  Wu-Dy  method  has  recently  been 
computed  by  Chaudhari,  Graczyk,  and  Charbnau(4)  for  the  Henderson  model.  Their 
results  confirms  our  original  conjecture  that  considerable  fine  structures  in 
this  function  will  remain  when  the  lattice  becomes  amorphous.  These  fine 
structures  are  due  to  coherent  scattering  of  clusters  of  atoms  with  short  range 
order.  In  order  to  complete  our  calculation  of  the  density  of  states,  we  need 

•  '  ’■  1  4 

in  addition  the  function  ,  which  is  the  fluctuation  in  the  energy  band  due 

to  disorder.  Since  computer  effort  in  this  direction  is  anticipated  to  be 

enormous,  we  are  now  trying  to  see  if  a  reasonable  analytic  approximation  is  • 

■  *’•••?  .  .  V 

possible.  ;< 

In  our  attempt  to  understand  the  electronic  structure,  we  have  also  investi*- 
gated  the  tight  -  binding  approach  to  the  problem.  We  have  calculated  the  density 
of  states  of  the  Henderson  modal  and  the  Polk  model  based  on  the  one  band  hamil- 
tonian  H  -  £  v|i)(i'J,  where  the  sum  is  over  all  nearest  neighbor  pairs  i  and  i'. 

In  order  to  calculate  the  density  of  states,  we  computed  the  first  ten  moments 
of  the  density  of  states  using  the  close-walk  counting  technique.  We  then  derive 
the  density  of  states  by  fitting  the  moments  with  an  expansion  in  Legendre  poly¬ 
nomials.  Since  the  Henderson  and  Polk  models  are  small  enough  for  the  hamilton<an 
to  be  diagonalized  numerically,  we  also  calculated  the  exact  eigen-values  and 
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eigen-states.  The  density  of  states  derived  in  this  way  agree  very  well  with 
those  derived  from  the  moment  method.  The  eigen-states  so  derived  exhibit  no 
localization.  In  order  to  produce  localization,  we  made  a  numerical  study  of 
the  Anderson  model  hamiltonian 

H  -  l  E|i)<i|  +  l  V|i><i' | 
i  ii' 

where  E  fluctuates  within  the  bounds  ±  W/2.  We  calculated  the  localization 


where  the  ai  are  amplitudes  for  some  particular  eigen-statejj.  In  a  similar  study, 
Edwards  and  Thornless  ^  show  that  for  perfect  diamond  lattice,  the  onset  of 
localization  is  less  sharply  defined  and  occurs  more  easily  than  in  the  Anderson 
theory.  For  the  disordered  Henderson  model  we  find  that  the  onset  of  localization 
is  also  not  sharply  defined  and  that  localization  is  the  same  on  the  average  to 
that  in  the  perfect  diamond  structure.  The  extension  to  more  realistic  models 
are  now  being  considered. 

(2)  Models  of  Steady  State  Electron  Injection Current  into  Dense  Media. 

We  have  completed  our  calculations  of  the  field  and  density  dependence  of  the 
steady  state  electron  current  injected  into  dense  argon,  hydrogen,  and  nitrogen 
gases.  The  results  are  included  in  a  preprint  attached. 
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Experiment 

Experimentally,  we  have  concentrated  on  continuing  our  experiments  on 
the  drift  of  electrons  in  dense  helium  gaa.  Our  results  are  now  in  accord  th 
previous  experiments,  the  precipitous  drop  we  reported  in  our  last  report  is  in 
fact  double  branching,  i.e.  a  high  mobility  branch  coexisting  with  a  low  mobility 
branch.  This  was  observed  before  in  hydrogen  but  never  in  helium  and  suggest 
a  very  long  life  time  for  bubble  formation  at  these  high  densities.  A  preprint 
covering  these  result  is  attached. 
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Abstract 

The  structure  of  the  bound  band  of  a  one-dimensional  disordered 
array  of  attractive  6-function  potentials  was  calculated  using  the 
method  of  Wu  and  Dy."*- 
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I .  Introduction 


In  a  previous  paper  (referred  to  as  I),  we  presented  a  method  for  cal¬ 
culating  the  electronic  density  of  states  in  disordered  systems.  We  showed 
that  the  method  can  be  easily  generalized  to  systems  with  structural  disorder. 

Here  we  applied  the  method  to  calculate  the  density  of  bound  states  for  an 
electron  in  a  one-dimensional  disordered  array  of  attractive  6-function  poten¬ 
tials  of  a  given  strength.  Other  methods  of  solving  this  particular  problem 
2“  6 

exist.  Our  purpose  here  is  to  illustrate  how  the  method  presented  in  I 
works  and  to  provide  insight  for  the  application  of  the  method  to  three-dimensional 
problems . 

Although  the  one-dimensional  disordered  array  is  not  a  topological  disordered 

system,  we  have  not  taken  the  advantage  of  the  sequential  numbering  of  the 

atomic  sites  to  reduce  the  problem  into  an  ordered  one.  We  also  do  not  use 

the  closely  connected  method  of  relating  the  energy  eigenvalue  to  the  number 

of  nodes  in  the  wavefunction.  Thus  our  method  of  calculation  can  be  readily 

extended  to  two  and  three  dimensions.  However,  as  a  basis  for  comparison,  we 

also  calculated  the  exact  density  of  states  by  the  node  counting  method  of 
o 

Lax  and  Phillips. 

We  shall  consider  only  arrays  having  short  range  order  as  defined  by 
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Gubanov  in  the  following  way.  The  distance  between  each  pair  of  neighboring 
atoms  is  taken  to  be  a(l  +  ey )  where  a  is  the  average  interatomic  spacing 
L/N,  L  being  the  length  of  the  chain  and  N  the  number  of  atoms  in  the  chain. 

The  factor  c  is  a  positive  number  less  than  one  (called  the  short  range  order 
parameter)  and  y  is  a  random  number  having  a  Gaussian  distribution  with  ^y)  =  0, 

=  1.  Makinson  and  Roberts^  have  computed  the  density  of  states  of  one¬ 
dimensional  arrays  of  6-function  potentials  with  c  ranging  from  0  to  01. 

The  results  presented,  however,  are  for  the  positive  energy  region  only.  In 
this  work  we  have  calculated  the  energy  level  distribution  in  the  negative 

energy  region  for  the  same  range  of  e. 
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II.  Method  of  Calculation 

We  use  the  Green's  function  formulation  presented  in  I  to  calculate  the  density 


(2.1) 


of  states.  The  Hamiltonian  we  consider  is  given  by 

H  =  2m  "  \  V(X  ' 

Jv 

where  VQ  is  a  positive  number  which  defines  the  strength  of  the  potential  and  t 
denotes  the  atomic  sites.  The  tight-binding  wavefunction  is  taken  to  be  the 

(2.2) 


eigenstate  l 

mV  mV- 

4>(x  -  .*)  =  (-sr)  exPt-  ^"lx  "  ll3 


*2  '  *2 


of 


h  *  b  -  V<*  •  •*> 


(2.3) 


with  eigenvalue 

2ft2 


(2.4) 


For  simplicity  we  set  w  equal  to  -1.  Using  (2.2)  one  can  easily  compute  the 
matrix  elements  of  the  total  Hamiltonian: 


mV, 


H.„  •  «S„,  +  2w  l  «p(-  -  t|  ♦  l.t"  -  Vl»- 

U'  U'  £,,^1 

Defining  the  position  l  by  a  dimensionless  parameter  t  =  l/a,  we  &et 

Hii'  =  w6u*  +  2w-  I®XP{"AH*'  "  •*!  +  ”  *'l^» 


(2.5) 


(2.6) 


l'7.t 


mV, 


where  A  =  _J>e.  Henceforth  we  shall  omit  the  '  sign  over  the  t's.  The  parameter 
ix2  , 

A  is  equivalent  to  e  in  Lax  and  Phillips'  notation. 

In  the  formulation  of  I  we  need  the  following  quantities: 

n  _  M-1  r  i(k-k').t  , 

Rkk,  =  H  6kk" 


V<k)  =  e 


ikL 


(2.7) 

(2.8> 


U(t,k)  5  JH  L  « 
L 


ikL 


-  V(k), 


(2.9) 


-  M-lr  ,T  .  i(k"-k’).l 


V.  =  N'  l  (I  +  RV'  e' 


•  U( £,k ' ) . 


(2.10) 


Typical  curves  of  R,  V,  U  and  W  are  shown  in  Figs.  1-3.  The  matrix  (I  +  R) 

needed  to  calculate  Wkk,  is  inverted  from  I  +  R  using  a  computer. 

Using  the  computed  values  of  V(k)  and  Wkk,,  we  calculated  the  Green's  function 

by  the  direct  summation  of  the  band  propagator  expansion: 

G(k,z)  =  l/(z-  V(k)  -  E(k,z)).  (2  11) 


where 


£<k,z>  =  iikk  +  Wkk,G(k',z)Hk,k  ♦ 


(2.12) 


We  truncated  the  series  for  the  proper  self-energy  E  at  the  second  term  as  shown 
in  (2.12)  and  solved  for  G  self-consistently  using  a  computer.  The  density  of 


states  was  calculated  from  G  by 


dN(E)/dE  =  n(E)  =  -  -  [  Im  G(k,z  -*•  E  +  i0+). 

IT  _ 


(2.13) 


III.  Results  and  Discussion 

In  Fig.  1,  the  curve  for  (R^J2  shows  considerable  structures.  These 

structures  are  closely  related  to  coherent  scattering  from  clusters  of  atoms 
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in  the  array.  P.  Chaudhari,  J.  F.  Graczyk  and  H.  P.  Charbau  recently  calculated 
the  quantity  lRlck«|2  th-ir  notation)  for  a  three-dimensional  random  network 

model  of  a  tetrahedrally  coordinated  amorphous  solid.  Their  result  shows  similar 
structures  in  the  | Rkk , | 2  curve  but  with  less  intensity  in  the  local  maxima  com¬ 
pared  with  our  one-dimensional  random  array. 

_  r  u  ikL 

In  Fig.  2,  the  V(k)  curve  (solid  line)  is  the  average  of  l  e  over 

L 

all  sites  i.e.,  it  is  the  energy  dispersion  of  the  electron  in  an  'average  crystal'. 
The  approximation  G(k,z)  «  l/(z-V(k))  gives  essentially  the  energy  band  of  the 
perfectly  ordered  array  and,  therefore,  is  a  poor  approximation  for  the  tail  states. 

The  structures  in  the  U(.t,k)  curve,  for  a  given  k  value,  presented  in  Fig.  3 
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represent  the  deviation  of  e*  from  V(k)  from  site  to  site  and  arise 

because  of  the  disorder.  Note  that  the  fluctuations  are  of  the  same  order  of 
magnitude  as  the  band-width  given  by  V(k)  alone.  In  Fig.  2  we  see  that  the 
transformation  of  U(t,,k)  to  (dashed  line)  reduces  the  fluctuations  to  values 

much  smaller  than  the  band  width  of  V(k).  This  induction  is  due  to  the  short- 
range  order  in  the  array  considered . 

In  Fig.  4(a)  and  (b)  we  show  (solid  lines)  the  calculated  density  of 
states,  dN(E)/dE,  for  a  chain  of  50  particles  with  short-range  order  parameter  e 
equal  to  0.05  and  0.1.  in  Fig.  5(a)  and  (b)  we  show  (solid  lines)  the  corresponding 
integrated  density  of  states,  N(E).  In  all  these  cases,  the  quantity  A  which 
determines  the  degree  of  overlap  of  the  neighboring  wavefunctions  is  taken  to  be 
10.  As  a  comparison,  we  plotted  in  the  same  figures  the  exact  results  obtained 
by  the  node  counting  method  of  Lax  and  Phillips  (histograms  in  Fig.  4  and  dots  in 
Fig.  5).  The  positions  of  the  band  edges  for  the  perfectly  ordered  system  are 
indicated  by  arrows  in  the  figures.  From  Figs.  4  and  5  one  can  see  the  gradual 
spreading  of  the  allowed  band  as  the  disorder  increases.  Because  of  the  small 
size  of  the  arrays  considered,  the  exact  density  of  states  shows  considerable 
fine  structures.  In  Fig.  6  we  show  the  exact  density  of  states  for  a  bigger 
sample  consisting  of  2000  particles.  It  is  evident  that  the  fine  structures 
diminish  as  the  sample  size  becomes  larger .  The  density  of  states  calculated 
by  our  method  using  the  50~particle  sample  is  smoothed  out  because  of  the  approxi¬ 
mation  used  in  summing  the  proper  self-energy  series .  The  finite  imaginary  part 
of  the  proper  self-energy  gives  rise  to  a  broad  Lorentzian  rather  than  a  delta 
function  at  each  eigenvalue  in  the  density  of  states. 

Ideally,  the  comparison  of  results  should  be  made  for  a  very  large  sample. 
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Because  the  computer  time  for  the  present  method  increases  with  sample  size  N 

3 

roughtly  as  N  ,  no  sample  larger  than  50  particles  were  attempted.  Nevertheless, 


Figs  4  and  5  show  that  on  the  average  the  band  tails  are  quite  well  reproduced 
by  our  present  calculations.  The  larger  discrepancies  actually  occur  at  the 
center  of  the  band.  The  reason  is  that  states  at  the  center  of  the  band  are 
quite  extended  (see  Fig.  7b),  and  as  the  disorder  increases,  the  overlap  of  the 
wavefunctions  can  become  quite  large,  thus  introducing  larger  errors  in  the 


tight-binding  approximation. 

We  conclude  that  using  the  present  method,  the  band  tails  of  disordered 
systems  with  short  range  order  can  be  calculated  quite  accurately.  However, 
because  of  the  large  computer  time  involved  the  application  to  three  dimensional 
systems  could  best  be  done  by  introducing  analytic  approximations  to  the  func¬ 


tionals  R^, 


and  W 


kk' 


for  an  infinite  system. 
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Figure  Captions 


Figure  1 


Figure  2 


Figure  3 


Figure  4 


Figure  5 


Figure 


Figure 


| R  , | 2  as  a  function  of  | k— k *  |for  50  6-function  potentials  with  short 
range  order  parameter  e  =  C.05. 

Typical  V(k)  curve  (solid  lir°)  and  the  real  part  of  W(k',k)  curve 
(dash  line)  with  k'  =  -4ir/5  for  a  linear  chain  of  50  6-function  poten¬ 
tials  with  e  =  0.05  and  overlap  parameter  A  =  10.  Th*  imaginary  part  of 

same 

W(k',k),  whi-u  is  not  shown  here,  has  the^order  of  magnitude  as  Re  W(k',k). 
Real  (solid  line)  and  imaginary  (dash  line)  parts  of  typical  U(.£,k)  curve 
with  k  =  -2tr/5  for  50  6-function  potentials  with  e  =  0.05  and  A  =  10. 
w  is  the  eigen-energy  of  the  bound  state  of  a  single  6-function  potential. 
Density  of  states  calculated  with  the  Dy-Wu  method  (solid  lines)  and  the 
node-counting  method  (histogram^)  for  50  6-function  potentials  with  A  =  10 
and  (a)  e  =  0.05  and  (b)  e  =  0.1.  Arrows  indicate  the  band  edges  for  the 
perfectly  ordered  system.  The  bound-state  energy  of  a  single  6 -function 
potential  is  chosen  to  be  zero. 

Integrated  density  of  states  calculated  with  the  Dy-Wu  method  (solid 
line)  and  the  node-counting  method  (dots)  for  50  6-function  potentials 
with  A  =  10  and  (a)  e  =  0.05  and  (b)  e  =  0.1.  Arrows  indicate  the  band 
edges  for  the  perfectly  ordered  system. 

Density  of  states  calculated  with  the  node-counting  method  for  2000 
6-*function  potentials  with  (a)  e  =  0.05  and  (b)  e  -  0.1. 

Eigenstates  for  (a)  E  =  4.66  x  10_4|w|,  (b)  E  =  -0.1  *  10_4|w|,  and 
(c)  E  =  -7.35  x  l0-4|w|  for  a  50-particle  system  with  e  =  0.1.  The 

* 

eigenstates  i|/r(x)  =  )'aa(E)+(x  “  .*)  are  obtained  by  finding  the  eigen- 
E  l  1 

vectors  of  the  Hamiltonian  matrix.  Only  the  amplitudes,  a^,  are  plotted 
in  the  figure. 
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Models  tor  Steady  Stata  Electron  Injection  Current  into  Denne  Had  la* 

P.  Snejtek,  N.  Silver,  I.  L.  Huang  end  K.  S.  Dy 
University  of  horth  Caroline  at  Chapel  Hill 
Chapel  Hill,  Korth  Caroline  27514 

At  at race 

Tvo  *9d«la,  the  Monte  Carlo  and  the  diffusion  nodal  for  calculating 
the  electron  current  due  to  injection  of  hot  electron#  into  dense  Mdla  are 
frtinttad.  In  the  diffusion  nodal,  the  current  can  he  derived  eiectly  fcy 
•jiving  the  continuity  of  current  equation  ounerlcelly  eatuning  a  fast 
energy  relaaetloo  tve- level  spate*.  Meuwor  a  very  elofle  analytic 
•agression  for  the  current  la  obtained  If  one  use#  a  strong  diffusion 
epprentneiloa  for  the  hot  electrons.  The  validity  and  errors  resulting 
free  the  apy foaloailon  are  dlecoseed.  When  the  dlffualoo  nodal  falls, 
u»  Ifepeee  a  elofltftad  Mania  Carlo  ends  l  for  calculating  the  current. 

Th*  nodal  In  a  •edification  of  on  earlier  oethod  hy  Toung  end  tredfcury 
(o  include  the  effect  ef  onlitple  ecattorlog  mi  the  influence  of  the 
inege  potential. 


10 


2 


1 .  Introduction 

Electron  injection  experiaenta  hove  played  a  powerful  role  in  determining 

12  3 

the  electronic  transport  properties  of  insulating  solids,  liquids  and  gases 
However,  up  until  recently,  there  la  atlll  a  lack  of  an  adequate  theoretical  model 
for  explaining  the  dependence  of  the  experimentally  measured  currant  on  the  applied 
electric  field,  the  characteristics  of  the  medium,  and  the  electron-medium  lnter- 
ictlon.  Perhaps  the  first  calculation  of  the  current-voltage  relation  resulting 
from  charge  tnjectfon  was  due  to  J.  J.  Thomson  4  .  Thomson  assumed  that  near  the 
cathode  the  Injected  electrons  have  a  random  velocity  v^  given  hy  thermal  equlllbtum 
with  the  gas  and  that  the  density  of  electrons  In  ths  gee  la  uniform  and  of  average 
value  n  *  This  Implies  that  there  Is  t  hack-diffusion  current  given  by  nev^/*  • 

In  the  presence  of  an  applied  field,  the  current  will  he  given  hy  J*nev  where  v 
la  the  drift  velocity.  Thus  If  JQ  is  the  total  injected  current,  v*  get  from 
Theme on*a  result  Jq-J***«q/*  •  **d 

J/V4S/<V*  #  (1) 

Theobald  '  end  Uefc  4  found  that  tg.  (1)  gives  reasonable  fit  to  experiment  only 
If  v  Is  reinterpreted  as  the  average  velocity  of  emission,  tekisrlan  7  Improved 
the  theory  farther  hy  taking  the  hnck-dlf feetem  current  te  he  a^/*  and  the 
measured  current  te  he  n^cv  dure  n^  end  naem  the  dens  Ulna  af  elactrona  near 
and  far  fren  the  cathode  respectively.  The  applicability  of  the  theory  presented 
above  1*  gelt*  restricted  even  In  Ite  modified  form.  Theobald  *  found  that  thera 
is  a  large  desetopency  between  theory  end  experiment  when  the  emits  toe  velocity 
be-* ernes  much  ls«s  them  the  drift  velocity.  In  such  s  case  hack-diffusion  loss 
de  teases  sad  Thomson *s  result  fella.  Another  situs  Hon  where  Ig.  (1)  may  he 
espected  to  fall  It  when  the  bech-dtffeelen  cone  la  oadsr  the  strong  Influence  of 
the  image  field  predated  by  the  Injected  electrons.  In  such  n  cese  *he  return  af 


20 


3 


electrons  Co  lh«  cathode  la  noc  given  by  back-diffusion  alone. 

Perhaps  Che  corrccc  approach  cc  Che  soludon  of  che  problem  would  be  Co  do 
a  Hence  Carlo  calculadon  including  all  cha  fields  experienced  by  che  elecCrons 
and  all  Che  laporcanc  acaccerlng  processes  encounter  by  Che  electron*.  This  Is  a 
very  dae  conaualng  calculadon.  One  alaplifled  approach  Co  chls  was  Bade  by 
Young  and  Iradbury  ®  who  calculeced  Cha  recurn  currsne  ansuaing  only  reflecdon 
of  elaccrona  In  chair  firac  encounter  with  gas  aeoas  or  aolccules.  They  found 
that  for  J/Jq<0.2  che  current  la  given  by 


where  \  is  che  as an  free  path,  E  is  che  applied  field  sod  is  cha  energy  of  cha 
injected  electrons.  Unforcunecely  che  validity  of  Eq.  (2)  is  difficult  to  access 
ee  chare  is  no  slap  Is  tray  of  estimating  che  effaces  of  auldple  acaccerlng. 

In  this  paper,  we  shell  present  two  theoretical  aodelt  for  che  sceedy  etece 
electron  injection  currents.  The  first  model  is  e  modification  or  tha  Young  and 
Iradbury  medal  to  account  for  the  effects  of  both  Image  field  and  auldple  acaccerlng. 
Th*  model  le  applicable  to  cases  where  the  electron  momentum  and  energy  relexedon 
rates  ere  alow.  The  second  model  le  a  two  fluid  model  in  which  the  electrons  are 
divided  into  e  hoc  end  e  thermal lied  component.  The  currents  are  derived  by  solution 
of  Che  continuity  ol  currents  equation.  This  model  is  sppllcable  when  tha  aomsntum 
and  energy  relaxation  races  of  che  electrons  are  rest. 

II.  Simplified  Monte  Carlo  Model 

In  this  see don  we  prsseac  s  simplified  Konte  Carlo  method  for  calculating 
che  current  1m  an  experimental  situation  depicted  in  Pig.  1.  An  electron  le  Injected 
Into  the  medium  with  energy  tQ.  In  the  taedUa  the  electron  experiences  the  total 
potential  ?(x)  consisting  of  the  image  potential  -  ^  and  the  applied  potential 
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-eEx.  Let  us  consider  first  the  case  ia  which  the  scattering  of  the  electron  is 
clastic.  In  our  simplified  modal  we  shall  male  the  following  basic  assumptions: 

(a)  The  scattering  of  an  electron  with  a  particle  in  the  medium  is  s-wave, 

(b)  The  direction  of  motion  of  the  electron  after  each  collision  is  inde¬ 
pendent  of  the  previous  collision. 

(c)  The  thermal  velocity  of  the  scatterers  in  :he  medium  are  negligible. 

(d)  The  trajectories  of  the  electron  between  collisions  are  linear  and 
their  path  lengths  a  distributed  according  to  exp(-e/A)  where  X  is  the 

momentum  exchange  mean  free  path. 

Assumption  (a)  can  be  easily  corrected  to  Include  other  partial  waves. 
Assumption  (b)  is  valid  if  the  average  collision  duration  of  an  electron  with 
the  scattering  center  is  much  shorter  than  the  time  between  collisions.  Assumption 
(c)  limits  our  presantatlon  to  mediums  with  extremely  low  temperature.  In  these 
cases  the  electron  loses  its  kinetic  energy  at  each  collision  only  through  the 
racoll  of  the  scattered  and  tha  amount  lost  is  simply  calculated.  This  assumption 
shortens  our  Monte  Carlo  calculation  considerably  but  it  is  by  no  means  necessary. 
At  higher  temperatures,  one  can  for  examply  assign  a  Haxvell-Bsltzmann  distribution 
of  velocities  to  the  particles  in  the  medium.  The  change  in  kinetic  energy  of  the 
electron  after  a  collision  can  then  be  computed  from  the  collision  kinematics.  In 
this  paper,  we  only  present  tha  T>0*K  results. 

Assumption  (d)  is  correct  if  the  mean  fres  path  of  the  electrons  is  so  short 
that  batwean  collisions  they  don't  pick  up  (or  lose)  enough  kinetic  energy  from 
the  field  to  deflect  their  trajectories  significantly.  In  this  case  the  problem 
reduces  to  a  straight  diffusion  problem  and  it  can  be  handled  simply  for  example 
by  solving  a  Boltzmann  aquation  treating  the  field  as  a  perturbation.  We  shall  be 
interested  in  casee  vhare  the  influence  of  the  image  and  driving  fields  are  strong 
so  that  assumption  (d)  has  to  be  corrected  to  take  into  account  trajectories  with 
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turning  points  and  the  fact  that  most  trajectories  will  be  bent  toward  the 

direction  of  the  field.  Since  such  calculations  would  require  a  tremendous 

amount  of  computer  time,  we  correct  for  the  effect  In  the  following  approximate 

way.  Let  us  consider  an  electron  located  at  a  position  x  to  the  right  of  the 

potential  maximum  x  as  shown  In  Fig.  1.  The  probability  that  the  electron  will 

M 

be  reflected  backward  by  the  fields  at  that  point  will  be  proportional  to  the 
8 

return  .one  , 

f  rV7^1/2, 

B<‘>  ’  }  (3) 

We  shall  define  a  reflection  coefficient,  r(x)  and  a  transmission  coefficient  t(x) 
for  an  electron  at  a  distance  x>x^  from  the  cathode  to  be  r(x)<*n(x)/4ir  and 
t(x)-l-r(x)  .  For  x^  ,  r(x)  and  t(x)  simply  exchange  roles.  In  our  model  we 
assume  as  In  (d)  that  the  trajectories  are  linear,  but  correct  It  approximately 
by  weighting  the  forward  and  backward  trajectories  according  to  t(x)  and  r(x) 
respectively.  Thus  our  simplified  Monte  Carlo  calculation  proceeds  as  follows: 

(1)  Electrons  are  injected  into  the  medium  with  energy  eq  and  with  an 

Initial  angular  distribution  f  (8  ,$  )  .  In  the  present  calculation  we  take 

0  0  u 

to  be  monoenergetic  and  f  (8^,^)  to  be  In  the  forward  direction  so  that  the  current 
can  be  correlated  with  injection  Into  vacuum  where  the  escape  cone  is  expected  to 
be  small. 

(2)  The  length  of  the  trajectory  s  of  a  given  electron  is  determined  by 
-Mnr ,  where  r  is  a  random  number  between  0  and  1,  and  X  is  the  momentum  exchange 
mean  free  path. 

(3)  At  the  end  of  the  trajectory  where  a  collision  event  has  occured,  the 
electron  losses  approximately  — (l-cos0)  of  its  kinetic  energy  to  the  recoil  of  the 
molecule,  m  and  M  are  the  masses  of  the  electron  and  the  scatterer  respectively 
and  9  is  the  angle  between  the  Incident  and  final  velocities  of  the  electron. 
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(4)  After  the  collision,  the  direction  of  motion  of  the  electron  with 
respect  to  the  x  axis  is  given  by  e.cos_1r  where  r  is  a  random  number  between  0 
and  1  .  The  choice  of  6  in  the  forward  or  backward  direction  is  weighted  according 

to  t(x)  and  r(x)  respectively. 

(5)  Taking  the  angle  c  f  step  (4)  as  the  new  initial  angle  we  repeat  the 
process  of  generating  the  tr Rectories  until  the  final  point  where  the  kinetic 
energy  becomes  less  than  the  potential  maximum  VM  .  Then  the  electron  is  considered 
transmitted  if  its  final  position  x  is  greater  than  Xjj  and  returned  otherwise. 

(6)  If  in  addition  there  are  inelastic  processes  in  the  medium,  we  need 

to  know  all  the  possible  modes  of  energy  transfer  and  their  respective  cross-sections. 
The  number  of  each  type  of  inelastic  events  is  then  generated  in  direct  proportion 

to  their  cross  sections. 

The  results  for  the  Monte  Carlo  calculation  of  the  current  as  a  function 

of  E/N  where  E  is  the  strength  of  the  applied  field  and  N  is  the  density  of  particle 

in  the  medium  is  shown  by  the  solid  curve  in  Fig.  2.  In  these  calculations  two 

atomic  masses  were  considered,  M-2  and  M-40,  and  only  elastic  scatterings  were 

—15  2 

Included.  The  scattering  aToss  section  was  taken  to  be  10  cm  and  the  density 

1021  atoms  per  cm3  .  The  results  are  independent  of  the  masses  to  within  the 

4 

statistical  fluctuation  due  to  the  finite  number  of  electrons  (10  )  injected. 

The  Young  and  Bradbury  result  is  shown  by  the  dashed  curve.  The  agreement  between, 
the  Young  end  Bradbury  result  and  the  Monte  Carlo  result  is  quite  good  at  high  field, 
for  the  density  and  cross-section  considered. Under  such  conditions  it  is  reasonable 
to  neglect  the  image  field  and  the  single  scattering  approximation  is  valid  because 
the  transmission  coefficient  rapidly  becomes  one.  At  lower  fields,  the  Monte-Carlo 

results  gives  lesser  end  lesser  currents  compared  with  the  Young  and  Bradbury  results. 
The  differences  arise  from  inclusion  of  the  image  field  and  multiple  scattering.  The 
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deviation  will  become  smaller  If  we  consider  finite  temperature  cases  as  part  of 
the  current  will  thermally  diffuse  across  the  potential  barrir 
III.  The  Diffusion  Model 

When  there  is  rapid  momentum  and  energy  relaxation  it  becomes  simpler  to 
calculuce  the  current  by  a  diffusion  model.  First  we  shall  discuss  the  injection 
of  monoenergetlc  electrons  into  the  medium  where  the  momentum  exchange  m.f.p.  is 
X  and  the  strength  of  the  inelastic  electron— medium  interaction  is  given  in  terms 
of  an  average  lifetime  t.  .  Because  we  use  only  a  single  relaxation  time 
the  system  can  be  characterized  by  a  two-components  fluid,  i.e.  hot  electrons, 
the  density  of  which  are  given  by  p^  ,  and  the  tharmalized  electrons  with 
specifying  their  density.  Because  7* (Jj+lt^"0  »  we  have  simply  'or  the  case  of 
planar  geometry: 


-  5L/i  )  - 

dTV 


Ph/r«0 


(A) 


where 


and 


VAtJ-  +  “h  <B*  ^2)eh 


(4a) 

(5) 


sher*  it’  ’  “tT&V'  (i,) 

In  these  equations  D  is  the  diffusion  constant,  y  is  the  mobility  and  « 
is  the  dielectric  constant  of  the  medium  which  we  shall  set  equal  to  1.  Since 
Eq.  (4)  is  decoupled  from  the  thermal  component  we  can  find  its  solution  inde¬ 
pendently.  One  boundary  condition  1b  given  by  the  current  balance  at  one 
scattering  m.f.p.  from  the  emitter.  One  part  of  the  electron  current  injected 
is  scattered  back  without  appreciable  energy  loss  and  is  given  by  the  Thomson's 
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term  p.  (A)<v  (A)>  .  <v  (.}.)>  Is  the  average  over  all  angles  of  the  x  component 
h  x  x 

of  the  velocity  at  A  in  the  presence  of  the  applied  and  Image  fields.  If  we  define 
<yx(A)>=c(A)yQ  where  ,  Is  the  average  emission  velocity,  then  eft)  must  approach 
the  Thomson  value  of  1/A  as  A  becomes  larger  or  of  the  order  of  x^  ,  the  position 
of  the  potential  maximum.  The  details  of  the  calculation  of  ciA)  are  shown  in 
appendix  1.  Another  part  of  the  current  at  A  is  that  diffusing  into  the  medium 

O* 

and  given  by  f  r"1p^(x)dg  .  Consequently  for  current  balance  at  A  ,  the  boundary 
A 

condition  is 

CD 

J0  -  c(A)phft)v0+  T_1jph(*)dx  (6) 

A 


where  Jq  is 
use  is 


the  current  supplied  by  the  emitter. 


jh(»> 


0  . 


Another  boundary  condition  we  shall 

(7) 


Integrating  Eq.  (A)  from  x  to  •  and  using  (7)  gives 

CD 

jh(x)  -  T“1|ph(x)dx  (8) 

X 

We  shall  show  that  the  measured  current  j  can  be  derived  from  Jh(x)  alone.  We 
start  by  discussing  the  T  -  0  K  case.  At  T  -  0  K,  the  diffusion  part  of  the  thermal 
current  given  by  Eq.  (5e)  /anishes  as  Dt-0.  At  x  -  x^  the  field  driven  current 

also  vanishes  as  the  field  is  sero  at  .  Thus  Jt(^)  ■  0  .  and  the  measured 
current 


3(T-0)  -  jh+jt  -  W  .  (9) 

For  finite  temperatures,  Jfc  is  nonvenlshing  at  x^  and  we  have  to  adopt  another 

method.  Since  J  «  J.  (x)  +  J..00  we  8et  «f*er  substituting  Eq.  (5a)  for  j  and  -  ~ 

lit  t  ax 

for  the  total  field  the  following  equation, 
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.  ,  ,  v  n  dpt  dV 

j-jh(x)  -  -Dt  S”  -Vt -3J  Pt 


Solving  this  equation  for  pfc  we  get 

x  j  (x1 )-j 

Pt(x)  -  exp-[(|it/Dt)V(x)3/  - exp[(yt/Dt)V(x')]dx'  +  p^O) 


(10) 


As  x-**°,  we  don't  expect  p^_  to  diverge,  but  since  V(x)-*—“,  we  must  require  that  the 
integral  in  Eq.  (9)  vanishes.  Thus  we  obtain 


J 


J.  (x)exp[eV(x)/kT]dx 

>  h _ 

/  exp[eV(x) /kT]dx 
0 


(ID 


where  we  have  made  use  of  the  Einstein  relation  p^/D^  “  •  We  shall  use  Eq.  (9) 

and  Eq.  (11)  to  calculate  the  current  and  as  mentioned  earlier,  they  depend  only 

on  Jh(x)  • 

We  now  discuss  the  derivation  of  P^(x)  from  which  can  be  calculated. 

Substituting  Eq.  (4a)  into  Eq.  (4),  we  get 


.2 

d  Pl 
D - £ 

h  di2 


-«h<E- 


dp 


4ex 


2'dx 


h  _(V 


2«x 


3  +  t)ph  '  0 


(12) 


Equation  (12)  can  be  solved  exactly  by  numerical  method.  Typical  solutions  are 

shown  in  Fig.  3.  Here  we  shall  present  a  simple  approximate  solution  and  discuss 

its  validity.  We  note  that  if  the  terms  in  Eq.  (12)  arising  from  the  image  field 

c^n  be  neglected,  then  the  remaining  equation  can  be  solved  analytically.  We  expect 

that  the  approximation  might  be  justified  for  x  above  a  certain  distance  x„  for 

dph  SD 

which  ®h~dx~  >>  **h  2  ^h  *  we  this  the  strong  diffusion  approximation 

4«X 

(SDA).  Estimates  of  xgD  are  given  in  detail  in  appendix  2.  The  SDA  solution  to 
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Eq.  (12)  is  8 Imply 

Ph(x)  -  Ph(X)  e"Y(x"X) 


(13) 


where 


.  r&K  x -w* 

Y  2Dh  ^2Dh;  X0  J 


and 


x0  -v 


i.- 


Applying  the  boundary  condition  (6),  and  -  -XvQ  „  we  find 


P.(A)  [C(X)+-^-]_1 

h  W 


(14) 

(15) 


(16) 


The  comparisons  of  the  SDA  solution  with  the  exact  solution  are  shown  in  Fig.  3. 
It  is  observed  that  the  deviation  ia  greater  for  electrons  of  lower  energy  and 
shorter  mean  free  path.  This  deviation  is  to  be  expected  since  under  these  con¬ 
ditions  the  diffusing  electrons  are  more  exposed  to  the  Influence  of  the  image 
field.  Their  random  velocity  becomes  also  more  sensitive  to  the  change  in  the 
potential  energy  and  the  electrons  undergo  more  scattering  events  c'.oaer  to  the 
emitter  where  the  retarding  field  is  high, 

We  can  now  calculate  the  SDA  expression  for  the  current  by  substituting 
(13)  into  (8)  and  then  use  (9)  or  (11).  We  shall  discuss  here  the  T*0K  case 
only  as  the  following  simple  analytic  expression  for  J  ollows  from  (9), 

4q  -  — jlC(A)  +  VY<VX> 

3yxq  3yxq  (17) 

A  typicel  current  versus  E/N  curve  is  shown  by  the  solid  curve  in  Fig.  4. 

For  illustrative  purposes  we  choose  the  medium  to  be  nitrogen  with  a  density  of 
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22  3 

10  molecules  /ca  .  We  assume  an  Inelastic  proceaa  In  which  an  electron 

loaea  0.25  «V  at  atch  colllalon  with  a  croaa  section  of  10**®  ce* .  and  v*  takr 

the  aoaintui  exchange  croaa  eectlon  <o  be  3a  10* 14  cm2.  The  electron*  an? 

Injected  with  energy  t  »laV.  The  field  dependence  of  the  current  enters  through 
1/2  ° 

•  (e/4«E)  .  gy  changing  the  applied  field  the  potential  maximum  can  be 

shifted,  In  thle  way  the  image  barrier  probaa  the  spatial  dletrlbutlon  of  the  hot 
electrone.  At  low  fielda,  x^  la  large  and  a  large  portion  of  the  hot  electron* 
relaxed  be. ore  reaching  the  potential  maximum.  Thle  gives  rise  to  an  exponential 
drop  In  the  current  given  by 


1 
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a,  a^O 


(A/3a0)/(O  ^-) 


(!•) 


At  large  fielda,  y  becomes  quite  email 


’  *  wW 


and 


1. 

J0 


me 

_ »- 

c(l)v0*v 


Cl*) 


Which  le  a  Imply  lq.  <l),  the  prediction  of  Itagbold  3  m*  Ueb4.  It  la  lnteraatlng 
to  Ice  chat  In  the  high  fielda  limit,  the  dlffuaion  modal  predicta  that  th# 
current  la  independent  of  the  loelaetlc  proceeeee.  The  reaeon  la  that  moat  electrons 
*r*  r*l**lng  beyond  the  potential  maximum  In  this  limit. 

We  have  elao  lnveetlgated  what  would  happen  if  the  diffusion  model  were  not 
applicable  and  the  Monte  Carlo  model  have  to  be  uaed.  The  results  are  presented 
In  Pig.  4.  in  the  high  field  region  (shown  by  dota)  agreement  In  ih-  order  of 
magnitude  of  the  current  can  be  achieved  only  by  using  9  .  Ajao"16^2 

9  m  1ft”l®_2  _  * 

1  »  •  Tk.  l.  th.  l..U«.c  lt  but  th, 

d*r’“4*”M  h**  ClW"  *««  "«  «•  «•»  lUi  b.h.,..,.  ,» 
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Aaauelng  the  therms Ilia t Ion  distinct  *o  •  100  A,  ws  find  the  critlcsl 

•  • 

distance  to  bs  17  A  snd  2f  A  for  lsV  and  0.3  eV  electrons  respectively. 

The  curves  In  Fl|.  3  shove  that  ths  SDA  solution  and  the  exact  solution  agree 
beyond  a  critlcsl  value  of  x  slightly  larger  than  that  set  lasted  here. 

Appendix  III. 

In  order  to  estimate  ths  error  in  the  alectron  scattering  m.f.p.  snd  the 

bot  electron  lifetime  derived  from  the  SDA,  we  make  the  following  gedankan 

experiment.  We  will  Inject  monoenargetic  alectron  currant  1  into  a  medium  at 

o 

T*OK  and  maasura  the  collected  current  which  is  determined  by  tha  value  of  the 

current  at  x^.  This  current  can  be  calculated  by  solving  numerically  Eq.  (12) 

for  assumed  values  of  the  actual  hot  electron  m.f.p.  JL-_  and  lifatlme 

ACT  ACT 

We  then  analyse  the  same  currant  using  ths  SDA  solution  (17),  and  find 

5DA 

tSDA*  The  results  are  shown  in  Pig.  6  where  two  cases,  (a)  c0  -  laV  and 
*  0*3  aV  are  shown.  Tha  m.f.p.  and  lifetime  assumed  are  labeled  X._ 

*0-  TACT  thoM  0btaln#<2  trom  SDA  Approximation  era  *SDA  and  tsda,  respectively. 
As  clearly  shown  the  SDA  underestimates  \  and  overestimates  t.  Tha  deviation 
is  larger  for  lower  Injection  energy  sad  smaller  m.f.p. 

Appendix  IV. 

In  terms  of  the  strong  diffusion  approximation  at  T  •  OK,  the  measured 
current  does  not  depend  explicitly  on  ths  energy  of  tha  Injected  electrons. 

If  the  therms 11 sat ion  distance  itself  was  energy  independent,  the  energy  distri¬ 
bution  of  electrons  at  tbs  barrier  maximum  would  be  identical  with  the  injected 
one.  It  would  then  bs  possible  to  correct  for  the  effect  of  the  energy  distri¬ 
bution  on  ths  error  in  m.f.p.  and  lifetime  by  weighting  each  with  ths  energy 
distribution  function.  However,  tha  assumption  of  equal  relaxation  for  electrons 
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of  different  energies  is  not  realistic.  It  is  rather  more  reasonable  to 
assign  equal  lifetime  to  all  electrons.  This  way  the  relaxation  distance 
becomes  energy  dependent  according  to  the  definition* 

*0  H  [<2£0/m)1/2ttT/3))1/2  (A10) 


In  the  constant  lifetime  approximation*  the  low  energy  electrons  have  shorter 
thermallsatlon  distance.  The  low  energy  electrons  are  more  strongly  attenuated, 
and  therefore,  their  contribution  to  the  measured  current  will  be  less.  Another 
energy  dependent  effect  comes  in  when  we  Introduce  the  image  potential.  Fig.  3 
shows  that  the  hot  electron  density  decrease  faster  in  the  exact  solution  than 
in  the  SDA  solution.  Thus  the  image  barrier  also  distorts  the  energy  distri¬ 
bution  of  current  at  x^.  The  two  effects  mentioned  above  were  Included  in  the 
estimate  of  errors  produced  by  using  the  strong  diffusion  approximation.  The 
error  was  estimated  from  the  following  procedure.  For  a  selected  value  of  the 
electron  m.f.p.  and  for  an  assumed  energy  distribution  of  injected  electrons, 
a  numerical  solution  of  the  exact  differential  equation  was  sought.  From  the 
solution  for  each  energy  a  contribution  to  the  hot  electron  current  is  calcu¬ 
lated  ualng  Bq.  (8).  The  total  current  is  then  found  by  averaging  over  all 
incident  energies  eQ  according  to  the  energy  distribution.  For  comparison  we 
also  analyzed  the  current  by  the  SDA  method  where  the  energy  distribution  of 
Injected  electrons  is  replaced  by  a  monoenergetic  one  with  sQ  equal  to  the 
average  energy  of  the  different  distribution  studied.  The  resulting  dis¬ 
crepancies  in  X  are  shown  in  Fig.  7a,  b,  and  c.  In  Fig.  7a  the  error 
derived  from  the  exact*  solution  and  the  SDA  are  shown  for  the  monoenergetic 
case.  The  result  will  serve  as  a  standard  for  evaluating  the  error  in  X 


whan  the  injection  energy  is  not  monoenergetic.  In  Fig.  7b,  the  assumed 

distribution  is  photoelectric  with  ^EL  m  /-/_  \2  ,  , 

dc  E  8ln(*e/Emax),  where 
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Emax  *  eav.ir/(iT^-4) . 
^  *  exp  (-e  /rAy)  • 


In  Fig.  7c,  the  distribution  Is  thermionic  with 


Figure  Cap t Ion t 

Fig.  1  Schematic  representation  of  electron  Injection  into  a  medium.  Due  to 

the  applied  electric  field  and  polarization  of  the  emitter  the  electron  potential 

a  1/2 

energy  has  a  maxipum  at  a  distance  x^  -  ■) 

Fig.  2  Comparison  between  results  of  Monte  Carlo  calculations  and  the  simplified 

theory  of  Young  and  Bradbury.  The  ratio  of  the  collected  to  t  injected  current 

j /.  as  a  function  of  the  ratio  of  the  applied  field  to  number  density  N.  The 
JO 

-15  2  21  -3 

parameters  are  6  ■  10  cm  and  N  -  10  cm  . 

Fig.  3  Distribution  of  hot  electrons  are  obtained  from  the  numerical  solution 
to  Eq.  12  (solid  line)  and  from  the  strong  diffusion  approximation  (broken  line) 

Q  °  «  « 

for  different  electron  m.f.p. :  1)5*,  2)10A,  3)20A,  4)40Af  and  for  the  energy 

of  injected  electrons  a)  lev;  b)0.3  eV. 

Fig.  4  Comparison  between  the  analytical  solution  obtained  from  the  strong 

diffusion  approximation  (solid  line)  and  results  of  Monto  Carlo  calculations. 

J/jQ  is  the  ratio  of  the  collected  to  the  injected  electron  current,  E/N  is 

the  ratio  of  the  applied  electric  field  to  number  density  of  the  medium.  Param- 

—16  2  —18  2 

eters:  strong  diffusion  approximation  0p  ■  3x10  cm  ,  oi  -  10  cm  ;  Monte  Carlo: 

triangles,  <jp  ■  1.9xl0-1^cm2,  -  2xl0_1®cm2.  circles,  0p  ■  4xlO*’1^cm2, 

—18  2  22  3 

<7^  "10  cm  .  Temperature  OK,  N  ■  10  molecules/cm  . 

Fig.  5  Backseat ter ing  coefficient  c(A)  of  hot  electrons  in  the  image  field 

as  calculated  from  Eq.  (A  4)  and  (A  5)  as  a  function  of  the  electron  m.f.p.  X  and 

energy  of  injected  electrons  e^. 
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Fig.  6  Evaluation  of  errors  produced  by  application  of  the  strong  diffusion 


approximation  to  the  analysis  of  the  measured  current  at  T«0K.  and 

tACT  are  rea*  of  the  hot  electron  m.f.p.  and  lifetime.  XSJJA  and  t£ 

are  the  corresponding  quantities  as  determined  from  the  strong  diffusion 


approximation  for  energy  of  injected  electrons  a)  cQ  -  leV,  b)  eQ  »  0.3  eV. 
Fig.  7  Evaluation  of  errors  produced  by  applying  SDA  method  to  the  analysis 


of  collected  current  when  the  distribution  of  injected  electrons  is  a)  mono- 
energetic  b)  photoelectric,  and  c)  thermionic. 
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Mobil  uy  of  Electrons  In  Dense  I'.eliua  Ga«* 

* 

J imi  A.  Jihnki  *nd  H.  Sllvir 
Department  M  Physics,  University  of  Kerch  Cerollne 
chapel  Mill,  Korth  Cerollne  27514 


tfe  iUvo  extended  the  rent*  of  electron  nobility  aeaaursaents 
In  dense  helium  |u  et  77.3  K  to  pressures  of  40  eta  end  et  160  K 
M  60  eta.  Coexisting  high  end  lea  nobility  branches  were  found 
et  both  temperatures,  corree ponding,  we  believe,  to  free  electron 
end  "bubble"  etetee. 

The  behavior  of  electron*  In  dense  hallux  gee  serves  *•  en  experimental 

1  2 

eodel  for  the  study  of  electrons  in  strvcturelly  disordered  eysteas.  '  A 
rapid  drop  in  electron  nobility  ss  s  function  of  density  was  suggested  by 
levias  and  Sanders'*  to  correspond  to  s  self-trapped  "bubble"  state. 

Recently  their  esperiaentel  werh  was  extended  by  Harrison  end  Springett*'* 
to  temperatures  above  the  heliua  critical  teaperature  (7  I  «  II  X)  shoving 

*  Research  supported  by  the  Advanced  Research  Projects  Agency  of  the  Depsrtaa 
of  Defense  and  aanitored  by  the  U.  S.  Army  Research  Office  -  Du rhea,  under 
Grant  Ko.  BA-A*DD-31-124-72-CiO,  end  by  the  Materials  Research  Center,  U.K.C. 
under  Cosinct  OAHC  15-67-C-0223  with  the  Advanced  Research  Projects  Agency. 


similar  rapid  decreases  in  nobility  as  a  function  of  fluid  density.  In 

contrast  to  the  self-trapped  model,  Eggerter  and  Cohen  interpreted  this 

data  in  terms  of  fluctuation  trapping  in  their  "pseudobubble"  model.  Due 

to  the  importance  of  distinguishing  between  the  forms  of  electron  trapping 

13  6 

in  these  theoretical  models,  *  *  it  was  felt  necessary  to  extend  the 

mobility  measurements  even  further  es  functions  of  temperature  and  pressure. 

Mobilities  arc  determined  from  the  measured  electron  drift  velocity  v  , 

d 

The  zero  field  electron  mobility  yu  is  defined  as  -  lim  v7E  where  I!  is 

^  0  E-*0  d 

the  electric  field  strength.  The  drift  velocity  was  measured,  in  the  present 

experiments,  by  use  of  a  double  gate  time-of-f light  technique.7  A  thin 

film  Al-Al^Oj-Au  diode  *  was  used  as  a  source  of  relatively  low  energy 

electrons  (1  eV  at  injection);  currents  of  about  10-13  to  10”11  amps  being 

successfully  swasured  in  the  present  apparatus.  Measurements  were  made  as 

a  function  of  pressure  at  77,3  K  and  160.1  K  by  immersing  the  experimental 

high  pressure  cell  in  dewars  of  liquid  nitrogen  and  frozen  ethanol.  Helium 

•  of  99.95Z  purity  was  passed  through  a  liquid-nitrogen-cooled  charcoal 

trap  and  pressurized.  Gas  of  99.9995Z  was  also  used  without  any  noticeable 

discrepancies  between  the  measured  drift  velocities.  Pressures  were 

determined  uelng  a  Helse  gauge  and  temperatures  were  monitored  with  a 

calibrated  copper-cons tantan  thermocouple.7  The  gas  number  densities  were 

calculated  from  an  empirical  equation  of  state, errors  estimated  as  being 

no  greater  than  +.02  x  1021  cm"3  at  77.3  K  end  +.04  x  1021  cm"3  at  160.1  K. 

Drift  velocity  measurements  were  taken  at  electric  field  strengths  of 

75  V/cm  to  200  V/cm  at  77.3  K,  but  only  at  100  V/c^  at  160.1  K  since  the  field 

dependent  data  at  the  lower  temperature  appeared  linear  in  this  range  for 

21  -3 

densities  greater  than  1.0  x  10  cm  (Fig.  1).  Each  point  on  the  drift 
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velocity  curves  is  the  average  of  over  12  experimental  traces  from  the 
spectrometer.  Over  six  harmonics  were  scanned  in  some  instances  for  a  given 
point,  whereas  no  less  than  two  harmonics  were  measured  for  any  given  point. 
The  precision  of  each  averaged  drift  velocity  is  about  5%  and  the  absolute 
accuracy  of  the  extrapolated  zero  field  mobilities  is  estimated  to  be  on 
the  order  of  10%. 

T'ne  zero  field  mobility  results  of  the  present  experiments  are  summarized 

in  Fig.  2  where  they  are  also  compared  with  the  recent  data  of  Harrison  and 

4  5 

Springett  for  both  helium  and  hydrogen.  *  As  can  be  seen,  our  data  at  77.3  K 

essentially  reproduces  the  earlier  data  over  the  range  at  which  it  was  taken. 

Also,  the  existence  of  two  mobility  branches  is  similar  to  that  previously 

observed  in  hydrogen.  The  mobility  branches,  we  observed,  coexist  over  a 

21  -3 

density  range  of  about  .7  x  10  cm  at  both  77.3  K  and  160.1  K.  The 
relative  amplitudes  of  the  drift  velocity  peaks  were  clearly  density 
dependent,  with  the  amplitudes  of  the  high  mobility  data  decreasing  with 
increasing  denoity  and  the  amplitues  of  the  low  mobility  branch  correspond¬ 
ingly  increasing.  The  amplitudes  of  each  branch,  at  constant  density,  were 
normalized  to  one  for  a  number  of  densities  and  are  plotted  in  Fig.  3.  It 
can  be  noted  from  the  electric  field  dependence  of  these  relative  amplitudes, 
that,  at  a  given  temperature,  the  high  mobility  branch  will  extend  to  higher 
densities  for  higher  field  strengths.  It  also  can  be  inferred  from  this 
data  that,  at  constant  field,  the  high  mobility  branch  will  extend  to  higher 
densities  a  the  temperature  increases. 

The  correspondence  between  our  helium  data  and  the  hydrogen  data  as 
shown  in  Fig.  2  (insert)  is  striking.  The  low  mobility  branch  in  hydrogen 
was  attributed  to  the  formation  of  ^  or  H  ions  by  the  high  energy  radio¬ 
active  source  used,  although  there  was  some  doubt  in  this  regard  by  the 
authors,  who  also  felt  that  self-trapped  "bubbles"  were  a  possibility."* 


In  che  present  experiment,  where  one  does  not  expect  the  formation  of 
negative  helium  ions11  since  a  low  energy  electron  source  was  used,  the 
low  mobility  branch  is  due  either  to  a  "bubble"  state  or  some  negative 
ion  impurity.  As  discussed  later,  we  rule  out  the  negative  ion  impurity 
as  a  possibility.  This  data,  then,  gives  support  to  the  contention  that 
the  low  mobility  branch  observed  in  hydrogen  is  also  due  to  bubbles. 

A  mming  then,  that  "bubble"  states  give  rise  to  the  low  mobility 
branch,  we  have  calculated  their  radii  at  these  higher  temperatures  using 
the  semi-classical  formula 

b  e  f1  .  9™ _ 1  CD 

p0  "  6nnR  4nR(*ifMItT)1/2 

where  n  is  the  viscosity,  R  the  bubble  radius,  n  the  number  density,  and  M 
the  reduced  mass  between  the  bubble  and  a  helium  atom  which  we  take  to  be  just 
the  mass  of  the  helium  atom.  We  obtain  a  radius  of  4.2  A  at  77.3  K  and  4.5  A 
at  160.1  K,  values  which  are  somewhat  too  low  when  one  considers  that  the  radii 
are  on  the  order  of  12  1  to  20  A  in  the  liquid.12*13  It  is  somewhat  doubtful, 
however,  whether  one  should  use  the  helium  mass  for  M  in  this  case  since  the 
"bubble"  mass  need  not  be  large  compared  to  that  of  the  helium  atom. 

For  several  reasons,  we  believe  that  impurities  do  not  play  a  measure- 
able  role  In  our  present  experiments.  In  the  first  instance,  we  used  helium 
gas  of  two  different  purity  levels  without  any  noticeable  differences  in 
drift  velocities.  In  the  erne  of  our  two  temperature  runs,  one  would  expect 
higher  impurity  levels  and  effects  at  160.1  K  than  at  77.3  K;  for  example,  an, 
oxygen  impurity  would  freete  out  at  77.3  K  but  not  at  160.1  K.  Such  an  impurity 
ion  would  then  give  a  low  mobility  signal  over  a  much  longer  density  range  than 
we  have  observed.7  Also,  in  regard  to  the  relative  peak  amplitudes  In  the  region 
„f  overlap  shown  in  Fig.  3.  the  signal  strength  due  to  a  negative  ton  tmpurttv 
would  increase  In  direct  proportion  to  the  density,  an  effect  which  has  not 
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been  observed. 


5 


in  order  Co  obtain  a  theoretical  estimate  of  the  zero  field  mobilities 

of  the  upper  branch,  we  have  used  a  modification  of  the  semi-classical 

14  5 

formula  for  gases  and  liquids  * 

c  1/2 

-  {4e/3a(2irmkT)  '  Xl+Bjn)  (2) 

where  m  is  the  electron  mass,  the  second  virial  coefficient  and  a  the 
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s-wave  scattering  cross-section.  At  densities  less  than  1.0  x  10  cm  our 

data  agrees  to  within  20%  of  these  theoretical  values,  but  for  higher 
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densities  this  agreement  becomes  gradually  poorer  until  at  3.0  x  10  cm 
these  theoretical  values  are  about  45%  too  high.  Trapping  due  to  fluctua¬ 
tions  (pseudobubbles)  would  lower  this  theoretical  estimate  to  essentially 
the  observed  value,  but  the  point  here  is  that  these  calculations  do 
suggest  that  we  are  dealing  with  free  electrons  in  the  high  mobility  branch. 

The  primary  implications  that  we  have  drawn  from  these  experiments 
are  that  we  have  observed  self-trapped  bubbles  in  dense  helium  gas  at  high 
temperatures  and  pressures  and  that  we  have  observed  the  coexistence  of 
free  electron  and  "bubble"  states  over  a  significant  density  range  (where 
quite  probably,  coexisting  free  electrons  and  pseudobubbles  lead  to  the 
high  mobility  branch  and  self-trapped  states  lead  to  the  low  mobility 
branch).  This  interpretation  further  suggests  that  we  are  dealing  with  a 
nonequilibrium  phenomenon  (because  of  the  branching)  implying  a  long 

relaxation  time  to  the  low  mobility  states.  This  is  of  some  interest  in 

8 

view  of  the  short  lifetimes  observed  in  the  liquid  and  in  view  of  the 
intermediate  mobilities  observed  at  lower  temperatures  where  fast  relaxa¬ 
tion  due  to  fluctuation  trapping  gives  rise  to  thermal  equilibrium:  between 
free  and  pscudobubble  states.  It  may  be  that  the  number  and  trapping 
cross  sections  of  the  states  observed  here  are  quite  strongly  dependent 
on  density  and  temperature,  giving  rise  to  a  longer  trapping  time  than 

52 


that  found  at  lower  temperatures.  If  one  were  to  obtain  a  longer  drift 
time  by  the  use  of  a  longer  drift  space  or  a  lower  electric  field,  one 
would  expect,  in  this  view,  that  the  free  electron  signal  strength  would 
decrease  with  a  corresponding  increase  in  the  "bubble"  signal.  This 
change  has  Indeed  been  observed  as  is  shown  in  Fig.  3. 

From  the  data  of  Harrison  and  Springett4  for  helium  at  18  K,  we  had 
expected  a  gradual  decrease  in  mobility  at  our  higher  temperatures  rather 
i-ian  branching.  It  is  possible  that  branching  also  occurs  at  these  lower 
temperatures  but  was  not  detected  by  the  single  gate  technique.  We  are 
presently  extending  our  experiments  to  lower  temperatures  in  order  to 
draw  some  consistency  between  these  two  sets  of  data. 

The  authors  wish  to  thank  J.  P.  Hernandez  for  many  helpful  discussions 
and  P.  Smejtek  and  T.  Carruthers  for  their  experimental  advice. 
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FIGURE  CAPTIONS 


Fig.  1 


Fig.  2 


Fig.  3 


Electron  Drift  velocity  versus  electric  field  strength  at  77.3  K 
in  a  pressure  range  of  10  Kg/an2  to  40  Kg/cm2.  The  drift  velocity 
is  nonlinear  at  100  V/cm  at  pressures  less  than  10  Kg/cm2. 
Zero-field  electron  mobilities  as  a  function  of  number  density  in 
dense  helium  gas.  The  dashed  lines  represent  the  earlier  data  of 
Harrison  and  Springett  for  helium.  The  inset  compares  our  77.3  K 
data  for  helium  with  the  hydrogen  dafa  of  Harrison  and  Springett 

at  31.7  K.  Note  the  existence  of  high  and  low  mobility  branches 
in  both  cases. 

Relative  signal  strengths  of  the  high  and  low  mobility  branches 
versus  number  density  at  77.3  K  and  160.1  K.  The  figure  for  77.3  K 
also  shows  the  effect  of  the  electric  field  on  the  relative  inten¬ 
sities.  The  data  for  160.1  K  was  taken  at  a  field  of  100  V/cm. 
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